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Abstract 

A Higgs particle produced in association with a Z boson and decaying into two 
photons is searched for in the data collected by the L3 experiment at LEP. All 
possible decay modes of the Z boson are investigated. No signal is observed in 447.5 
pb -1 of data recorded at centre-of-mass energies from 192 GeV up to 209 GeV. 
Limits on the branching fraction of the Higgs boson decay into two photons as a 
function of the Higgs mass are derived. A lower limit on the mass of a fermiophobic 
Higgs boson is set at 105.4 GeV at 95% confidence level. 
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Introduction 



The Standard Model of the electroweak interactions allows the decay of a Higgs boson h into 
a photon pair only at the one loop level. The branching fraction of this decay is small pQ. 
It amounts to about 0.1% for Higgs masses m h between 80 and 110 GeV. Several extensions 
of the Standard Model predict enhancements of this branching fraction For instance, an 
appropriate choice of the parameters of the Two Higgs Doublet Models of Type I [Hj, predicts 
the lightest CP even Higgs not to couple to fermions at tree level. Such a Higgs is expected 
to decay dominantly into a pair of photons if its mass is below 90 GeV In Higgs triplet 
models, one of the neutral scalars may have a large branching fraction BR(h — * 77) and be 
produced at LEP with rates comparable to the Standard Model ones [5 3 . At higher masses, the 
branching fraction BR(h — > 77) decreases to the advantage of the W+W+* mode. 

This letter presents the search for a Higgs boson produced in association with a Z boson 
through the process e+e - — > Zh, followed by the decay h — > 77. All decay modes of the Z 
boson are investigated. 

The high energy data sample collected by the L3 detector [2] at centre-of-mass \fs up to 
209 GeV is analysed. Results from partial samples were previously reported by L3 0|B] and 
other LEP Collaborations jH]. 

Data and Monte Carlo Samples 

We analyse data collected with the L3 detector during the years 1999 and 2000 at centre-of-mass 
energies y/s = 192 — 209 GeV, for a total integrated luminosity of 447.5 pb _1 . The data are 
grouped into ten samples whose average centre-of-mass energies and corresponding integrated 
luminosities are listed in Tabled Results from 176 pb _1 of integrated luminosity collected at 
y/s = 189 GeV are included in the final results. 

The Standard Model Higgs production cross section is calculated using the HZHA gen- 
erator PU]- Signal Monte Carlo samples are generated using PYTHIA [II] for Higgs masses 
between 50 and 120 GeV. These samples comprise between 500 and 2000 signal events depend- 
ing on the search channel. For background studies the following Monte Carlo programs are 
used: KK2f [12] (e+e" -»■ qqfr), e+e" -> // + /x~(7), e+e" -»■ r+r"^)), PYTHIA (e+e" -> ZZ 
and e+e" -> Ze+e"), KORALW [13] (e+e" -> W+W"), PHOJET [ft] (e+e" -> e+e'qq), KO- 
RALZ [IS] (e+e" -»■ ^(7)), GGG jTE] (e+e" -> 77(7)), BHWIDE [II] (e+e" -> e+e- (7)), 
TEEGG [U] (e+e" -> e+e"7), DIAG36 [U] (e+e" -> e+e'e+e") and EXCALIBUR [20] for 
other four-fermion final states. The number of simulated events for the most important back- 
ground channels is at least 100 times higher than the corresponding number of expected events 
in the data. 

The L3 detector response is simulated using the GEANT [21] program, which takes into 
account the effects of energy loss, multiple scattering and showering in the detector. The 
GHEISHA [22J package is used to simulate hadronic interactions in the detector. Time depen- 
dent inefficiencies, as monitored during the data acquisition period, are also simulated. 

Analysis Procedure 

The analysis aims to select events with isolated photons and a Z boson. Three final states 
are then investigated, according to the Z boson decay. They are denoted as qq77, ^z/77 and 
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£ + £~77) with I = e, /i, r. The selections for each final state are described in the following and 
proceed from a common photon identification. 

Photons are selected from clusters in the BGO electromagnetic calorimeter with an energy 
greater than 1 GeV and a shower profile compatible with that expected for an electromagnetic 
particle. The ratio of the energies measured in a 3x3 and a 5x5 matrix, centred on the most 
energetic crystal, must exceed 0.95. The energy deposit in the hadron calorimeter is required 
to be below 20% of the energy measured in the electromagnetic calorimeter. 

No track should point to the clusters within 50 mrad in the plane perpendicular to the 
beam axis. Only clusters in the polar angle ranges 25° < 9 < 35°, 45° < 9 < 135° and 
145° < 9 < 155° are considered, well within the coverage of the barrel and end-cap regions of 
the BGO electromagnetic calorimeter. The distribution of the polar angle of the most energetic 
photon is shown in Figure ^ 

These analyses require at least two photons. Photons from the decay of a heavy resonance 
are expected to be relatively energetic, hence the energy of the hardest photon has to exceed 
10 GeV and the energy of the second most energetic photon has to exceed 6 GeV. 

The qq77 Final State 

Candidates in the qqy7 final state are characterised by a pair of isolated photons accompanied 
by two jets. The qqy7 selection starts from a preselection of hadronic events: high multiplicity 
events are required to be balanced and their visible energy E vis /^fs has to be larger than 0.5. 
The background from two-photon interaction events is reduced requiring the energy in a 30° 
cone around the beam pipe to be less than half of the visible energy. The preselection yield 
is reported in Table El Figure El shows the comparison between data and the Monte Carlo 
expectations for the distribution of E vis /^/~s for the preselected events. 

Events which contain at least two photons are selected. In addition, photons coming from 
neutral hadron decays are rejected by requiring the energy in a 10° cone around the photon 
direction to be less than 2.5 GeV, and that in a 20° cone to be less than 4.5 GeV. The number 
of charged tracks and calorimeter clusters in a 20° cone around the photon direction must be 
below four. The opening angle between the photons must be larger than 50°. All other particles 
are clustered into two jets using the DURHAM jet algorithm [23] and no jet is allowed within 
25° around any photon. 

The energy spectrum of the most energetic photon, normalised to the beam energy, before 
the application of any selection requirement on the photon energies, is shown in Figure Ok- 
Figure 0b presents the distribution of the recoil mass against the di-photon system after the 
cuts on the photon energies. 

Finally, the recoil mass against the di-photon system is required to be consistent with the 
Z mass within 15 GeV. The efficiency is 40% for a Higgs boson mass of 100 GeV produced at 
v/i = 192 GeV, and 47% for a Higgs of 110 GeV produced at ^ = 208 GeV. The selection 
yield is presented in Table El 28 events are observed in the data, with 31 expected from Monte 
Carlo, mainly from the e + e~ — > qq(7) process. 

The event with the highest value of the di-photon invariant mass is displayed in Figure 01 
It was collected at a/s = 205.1 GeV and its di-photon invariant mass is 111.8 ± 1.0 GeV while 
the recoil mass against the di-photon system is 87.1 ± 0.8 GeV. 
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The 1/1/77 Final State 

The signature of the z/z/77 final state consists in two photons and missing energy. Events 
are selected that have an identified photon pair, no charged tracks and an additional energy 
below 10 GeV. To ensure that the missing momentum is well contained in the detector, the 
absolute value of the cosine of its polar angle must not exceed 0.96. To reduce contributions 
from the e + e~ — > 77(7) process and from double radiative events with final state particles 
escaping detection, the photon acoplanarity is required to exceed 3°. The distribution of this 
acoplanarity for the data and the Monte Carlo predictions is presented in Figure OJl The total 
transverse momentum of the di-photon system must be greater than 2 GeV. 

Figure shows the distribution of the recoil mass against the two most energetic photons 
after selection requirements on all the other variables. As final selection criterium, this mass 
has to be consistent with the Z boson mass within 15 GeV. The efficiency is 47% for a Higgs 
boson mass of 100 GeV produced at yfs = 192 GeV and 51% for a Higgs of 110 GeV produced 
at y/s = 206 GeV. The number of selected events in data is 9 and 9.2 events are expected from 
the e + e~ — ► vvi^j) process. Other backgrounds are negligible. 

The £ + £~^ Final State 

The £ + £~77 final state has the characteristic signature of a photon pair and a lepton pair. Its 
selection proceeds from low multiplicity events with two identified photons and an associated 
lepton pair, selected as follows. 

Electrons are identified from clusters in the electromagnetic calorimeter with an energy 
greater than 3 GeV and associated to a charged track. The energy deposit in the hadron 
calorimeter must be consistent with the tail of an electromagnetic shower. Less than 3 GeV 
are allowed in the electromagnetic calorimeter in a 10° cone around the electron direction. To 
increase efficiency, events with just one identified electron are also accepted. 

Muons are identified from tracks in the muon chambers with a distance to the interaction 
vertex in the r — plane below 300 mm and a momentum above 3 GeV. The calorimetric 
energy in a 10° cone around the muon direction must not exceed 3 GeV. Events with one muon 
and one minimum ionising particle in the calorimeters are also accepted, as well as events with 
a single muon. The background from cosmic rays is rejected by requiring at least one hit in the 
scintillation counters in a 5 ns window around the beam crossing time. 

Taus are identified as jets with one or three tracks in a 10° cone with an energy above 
3 GeV. The energy in the 10° — 30° cone around the tau direction has to be below 30% of the 
energy in the 0° — 10° cone. Events with only one identified tau lepton are also accepted. 

The preselection of ^ + ^~77 events yields the results listed in Table H3 Double radiative 
di-lepton events are the dominant background and are rejected further by requiring the energy 
of the most energetic lepton to be below 80 GeV. 

The distribution of the energy of the second most energetic photon normalised to the beam 
energy is presented in Figure |Hk for the preselected events. Figure 0) shows the recoil mass 
against the photons after all selection criteria but the recoil mass. 

The final selection requirement imposes the recoil mass to be consistent with the Z mass 
within 15 GeV. At centre-of-mass energies below y/s = 202 GeV the presence of two identified 
leptons is enforced. Their invariant mass is required to be between 81 and 101 GeV and the 
selection criterium on the recoil mass is relaxed. 

The efficiency varies from 31% for a Higgs boson mass of 100 GeV produced at y/s = 
192 GeV, to 43% for a Higgs of 110 GeV produced at y/s = 208 GeV. The yield of this 
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selection is presented in Table El 7 events are observed in the data, with 8.0 expected from 
Monte Carlo, mainly from the e + e _ — ► £ + £~ (7) processes. 

Results 

No significant excess indicating the production of a Higgs boson decaying into two photons 
is observed in the data. The confidence level [21] for the absence of a Higgs signal is then 
calculated from the reconstructed di-photon invariant mass as the final discriminant variable. 
This distribution is shown in Figure [7^, for the data analysed in this paper, collected at yfs = 
192 — 209 GeV and in Figure [7)d including the data collected at */s = 189 GeV. 

The calculation of the limits takes into account systematic uncertainties of 2% on the signal 
expectations and 8% on the background. The signal uncertainty follows from the Monte Carlo 
statistics that also accounts for a 4% uncertainty on the background. Another uncertainty of 
7% is assigned to the background normalisation for hadronic events with photons. A variation 
of 2% of the calorimetric energy scale has little effect on the limits. The effects of the energy and 
angular resolutions of the photons and the systematic uncertainty on the integrated luminosity 
are also found to be negligible. 

Figure |H1 presents the upper limit on the branching fraction BR(h — > 77) as a function 
of the Higgs mass, assuming the Standard Model cross section for the Zh production. The 
expected limit is also shown together with the theoretical prediction for a fermiophobic Higgs 
boson as calculated with the HDECAY program j2Hj. Previous L3 results [Hj are included in 
the calculation of this limit. The observed limit for BR(h — > 77) = 1 is 114 GeV. 

The lower limit on the mass of a fermiophobic Higgs boson is set at 

m h > 105.4 GeV at 95% confidence level, 
to be compared with the expected mass limit of 105.3 GeV. 
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1Q1 6 195 5 199 5 291 7 203 8 205 1 206 3 206 6 208 208 6 


Luminosity (pb _1 ) 


29.4 83.7 82.8 37.0 7.6 68.1 66.9 63.7 8.2 0.1 



Table 1: The average centre-of-mass energies and corresponding integrated luminosities. 





N D 


N B 


qq(7) 


WW 


Ze+e" 


ZZ 


Preselection 


17719 


17739.2 


11364.2 


5876.7 


139.4 


358.9 


Selection 


28 


31.0 


30.5 


0.2 


0.1 


0.2 



Table 2: Number of events, iV D , observed in data by the qq77 selection, compared with the 
Standard Model expectations, N B - The Monte Carlo breakdown in different processes is given. 





N D 


iV B 


e + e (7) 


fi+fi (7) 


T + T (7) 


4 fermion 


Preselection 


738 


751.3 


541.7 


46.2 


50.4 


113.0 


Selection 


7 


8.0 


4.2 


1.8 


2.0 


0.0 



Table 3: Number of events, iV D , observed in data by the £ + £ 77 selection, compared with the 
Standard Model expectations, N B - The Monte Carlo breakdown in different processes is given. 
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Figure 1: Distribution of the polar angle, 7l of the most energetic photon for data and back- 
ground. A Higgs boson signal with mass = 105 GeV is superimposed with arbitrary nor- 
malisation. All Z final states are combined. The selected regions are indicated by the arrows. 
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Figure 2: Distribution of E vis /y/s after the hadronic preselection for data and background for 
the qq77 final state. A Higgs boson signal with mass = 105 GeV is superimposed with 
arbitrary normalisation. The arrow indicates the value of the cut. 
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Figure 3: Distributions for the qq77 final state of a) the energy of the most energetic photon 
normalised to the beam energy and b) the recoil mass against the di-photon system in data, 
background and for a 105 GeV Higgs boson signal with arbitrary normalisation. The arrows 
indicate the values of the applied cuts. 
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Run # 852001 Event # 2399 




Figure 4: The e + e ^Zh^ qq77 candidate with the highest di-photon invariant mass. 
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Figure 5: Distributions for the final state of a) the acoplanarity of the 77 system and b) 
the recoil mass against the two photons in data, background and for a Higgs boson signal with 
mass mh = 105 GeV with arbitrary normalisation. The arrows indicate the values of the cuts. 
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Figure 6: Distributions for the £ + £~ , yy final state of a) the energy of the second most energetic 
photon normalised to the beam energy, for the preselected events, and b) the recoil mass against 
the two photons. Data, background and a Higgs boson signal with mass = 105 GeV and 
arbitrary normalisation are shown. The arrows indicate the values of the cuts. 
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Figure 7: Distributions of the reconstructed di-photon invariant mass for all final states com- 
bined, after the final selection. Data at a) y/s = 192 - 209 GeV, and b) y/s = 189 - 209 GeV 
are shown together with the background and a Higgs boson signal with mass = 105 GeV. 
The Standard Model cross section and a BR(h— > 77) = 1 are used. 
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Figure 8: Excluded values at 95% confidence level of BR(h — > 77) as a function of the Higgs 
mass, in the assumption of the Standard Model production cross section. The expected 95% 
confidence level limit and the theoretical prediction are also presented. 
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